function y = func_nxy_from_nz(nz, CONSTS)

    eps_a = CONSTS.eps_a;

    y = sqrt_from_negative_num(eps_a-nz.^2); 
    
end